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The physical origin of pairing in high temperature superconductors has not been de- 
termined as yet. Therefore, it is very important to formulate a language describing those 
systems, which is independent of the detailed underlying electronic structure, but still con- 
crete enough to allow for an extensive experimental test. We believe that the model proposed 
by Anderson and his coworkers [1] and based on the pairing enhancement due to an interlayer 
Cooper-pair tunneling provides such a crucial feature, on the basis of which one can con- 
struct an appropriate quantitative language. Universal features of this model not discussed 
so far are detailed below and compared with experiment. 

In this paper we derive analytically the space profile of the superconducting gap param- 
eter assuming only that we have a k-dependent pairing potential Vkk' in the two spatial 
dimensions (for a square planar configuration), and a coherent (Josephson) pair tunneling 
between the planes [1,2]. Moreover, we analyze explicitly the periodic arrangement of groups 
of planes, each of them containing p tightly spaced identical Cu02 planes, and separated 
by a larger distance, but in contact with each other also by, albeit weaker, Cooper pair 
tunneling. Such a configuration is drawn schematically in Fig. la for the high temperature 
superconductor. The system is thus characterized by the intragroup (T^(k)) and intergroup 
(T(k)) tunneling amplitudes. Hence, our approach can be applied to both mono- and multi- 
layer high-T c systems and thus generalizes previous treatments [1,3,4]. The calculated below 
space profile of the superconducting gap is shown schematically in Fig. lb. Our results for 
the transition temperature T c for either Fermi, or Luttinger or statistical spin liquid, the 
three universality classes appearing for the interacting fermion systems [5]. 

We label the single electron states by the in-plane component of the wave vector k = 
(k x , k y ) and by spin quantum number a = ±1, as well as by the two discrete numbers (s,j), 
where j — 1, . . . ,p characterizes the plane position within the group of tightly spaced planes, 
and s = 1, . . . , N labels the groups. The simplest Hamiltonian describing such a nontrivial 
structure can be written as follows: 
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The first two terms describe the BCS-type Hamiltonian for iV • p independent planes, the 
second and the third express the pair tunneling (with the matrix element T^(k)) between the 
planes in the same group, separated by (I — 1) superconducting planes in between, whereas 
the last term represents the tunneling between neighboring groups. We have assumed that 
the intergroup tunneling matrix element T(k) is independent of the pair of planes selected 
in those two groups, since the detailed balance within the group must be achieved much 
faster than between the groups. 

In the following we make the mean field (BCS-like) approximation, i.e., assume that 
the superconducting gap is characterized by the nonvanishing average dj^ =< ^kt^-kj >> 
which will not depend on group index s, as each group is in identical environment and 
N — > oo (this, however, is not the case for the individual planes within the group!) In effect, 
the gap function can be defined as follows 

p i i+'^p i<p 

A ilk = E2T(k)d ilk + -Evkkfd i> k'+ E r»(k)d i+l>k + "E" ^(kK-i*. (2) 
j=l iv k' 1=1 1=1 

The selfconsistency condition for the mean-field solution requires that 

dj,v = AjMXj, k- (3) 

The generalized susceptibility XjM i n this universal relation depends on the nature of the 
normal state. In the case of the Fermi liquid (FL) it takes the usual form [3]. 

X,k = ^-tanhf^, (4) 



where the quasiparticle energy is £j k = y (ek — p) 2 + | Aj ik | 2 . In the case of the spin-charge 
separated Luttinger liquid (LL), i.e., if we have the linear dispersion relation e^— p ~ (k—kp) 



for bare particles and two Fermi velocities corresponding to charge (v c ) and spin (v s ) degrees 
of freedom, the function Xj,k is [4] 

where E jM = \[v+ + yjvi + 4| A jik ] 4 ], h c , a = ^\x 2 - (v c>s k) 2 \, and v± = (v 2 ± v 2 )k 2 . Fi- 
nally, for the strongly correlated liquids with non-standard anticommutation relations for 
the fermion operators [5] the mean-field approximation of the BCS-type leads also to Eq. 
(2), in which rhs is renormalized only by the factor A k = 1 — n k /2, which near the Fermi 
energy has the value |. 

We take into account only the coupling between the nearest neighbors, i.e., assume that 
TW(k) = T(k)5 M± i. This approximation means that we can seek a solution in a simple 
exponential or trigonometric form. In that situation the gap A jk is determined by the 
system of universal difference equations 

1 p 
A ik - T7 XVk k /Xi,k'Ai, k / = T ( k )X2,kA 2 , k + 2T(k) ^ Xj',kAjv !k , 
iy k' j'=i 

1 p 
A i,k - -T7 Vkk'Xj,k'Aj ik / = T(k)(x i+ i k A i+lik + X j-i,kAj-i,k) + 2T ( k ) H Xj',kA /ik , 
iV k' j'=i 

1 p 
A p , k - — VWxp.k'Ap^ = T(k)xp-i, k A p _i ik + 2T(k) ^ xy >k Aj/ ik . (6) 

iV k' j'=i 
The shape of the spatial profile of A J;k (i.e., its dependence of j) should not depend on 

the system temperature, but only on the form of the boundary conditions. Therefore, we 

solve this system of equations close to the transition temperature T c , where the generalized 

susceptibility Xj,k can be linearized in the following way 



tanh (f £k) 



Xj^ 2e V J ' V +0(A£ k ), (7) 
for FL, and 

for LL. To zero order in the gap parameter the generalized susceptibilities do not depend 
on the layer index j. With this expansion, the system (6) represents a system of linear 



and inhomogeneous difference equations with respect to index j. In order to solve them we 
introduce two additional fictitious layers labeled by j = and j — p + 1 [6]. The natural 
boundary conditions are expressed in terms of the vanishing gap on these limiting layers, 
i.e., we set 

A 0>k = A p+lik = 0. (9) 

These boundary conditions express the physical equivalence of all p layers within the group, 
i.e., we do not introduce the surface layers, which differ from those inside. In such a situation, 
the system of equations (6) can be rewritten in a more symmetric form, both for the FL 
and the LL cases, namely 

- T7 E WkA,k' = T(k)xk(A j+lik + A 3 -_i jk ) + Z k , (10) 
iV k' 

where Z k = 2T(k)xk E/=i Af , k , and j = 1 . . .p. 

We solve first the homogeneous part of Eq. (10) (i.e., put Z k = 0). The spatial de- 
pendence of the solution for the gap can be taken in the form A J)k = A k ± ^e ±M *- J , where the 
constants A k ^ fulfill the set of self-consistent equations 

A k ±} (l -2T(k) Xk cos«) = ^E^k'Xk'AW (11) 

iV k' 

Thus, the general solution of the homogeneous part is 

A jik = A k +) e^' + A^e- iaj , (12) 

where the wave vector of the oscillating gap is determined from the boundary conditions 
(1), and takes the form 



( i i \ 



e ia(p+l) e ~ia(p+l) 



0. (13) 
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The vanishing determinant of the above matrix provides a nontrivial solution only when 
a = nn/(p + 1), where n is an integer. Substituting back the value of a to (13) we find 
A k + -* = — = A k ; so the solution (12) takes the explicit form 



A j)k = 2iA k sin^-^j. (14) 

Physically, the admissible solution is that with n — 1, as the solution with no nodes inside 
p layers should have the lowest energy. The shape of the spatial dependence of the gap has 
been drawn already in Fig. lb. Finally, the gap magnitude A k is determined from Eq. (10) 
with Z k = 0, which reduces to the integral equation: 

A k = — r- — V Vkk'Xk'A k < (15) 

l-2T(k) Xk cos(^ T )iV^ 

The solution of the full Eq. (10) is found by superposing the general solution(14) of the 
homogeneous part of (11) with a particular solution of the full equation, i.e., by taking 

A,, k = A k + A® sin , (16) 

where A k ^ = 2iA k , and A k is treated as a small perturbation. The part A k does not 
depend on j because the periodic boundary conditions can be taken for iV groups of layers. 
Substituting the solution (16) into (10) and omitting higher-order terms (i.e., those ~ A k • 
T(k)), we have the equation for A k ^ 

k l-2[T(k) / (p) + T(k)cos(^ T )]xk^& k ' 
The equation for A k is 

^■d^s^ (18) 

where f(p) = sin(7rp/2(p + 1))/ sin(7r/2(p + 1)). 

The physical T c is provided by the larger transition temperatures determined from Eqs. 
(17) and (18). Hence, it is always determined from Eq. (17). This feature of the theory 
illustrates the importance of the intergroup hopping T(k) and formally appends the ideas of 
[1]. Also, it reduces to that of Chakravarty et al. [1] in the limit of f = and p — 1. Most 
importantly, our result is the same for both Fermi and Luttinger liquids, as well as for the 
liquids with a generalized exclusion statistics described in [5,7]. 
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To discuss the implications of our approach let us consider first the monolayer case 
(p = 1), i.e., a material with the interlayer hopping T(k) = (t') 2 /t, for which t is the 
magnitude of the bare single-particle hopping in the plane. In that case a simple estimate 
of the denominator of Eq. (17) is: 

l-f(k) Xk ~l-/--^-, (19) 

l K B 1 C 

where r is a numerical constant of the order of unity. Hence, the tunneling can strongly 
enhance the pairing potential, since ksT c ~ r(t') 2 /t [1]. 

To calculate explicitly the value of T c we consider explicitly the p > 1 case with a 
constant pairing potential Vkk' = V in the energy regime of the width uj d <C /i. Introducing 
the density of states N(0) at the Fermi level we find from (17) 

1 = ™<°> 12 * i-w<pn£*)mmvt - m 

where a is the ratio of intergroup to intragroup tunneling amplitudes. In the extreme 
situation, when the interplanar tunneling provides a dominant contribution to the effective 
pairing, we can make the assumption ud/IcbTc <C 1 and thus approximate tanh(rr) pa x. 
Here, we provide the T c derivation for the Luttinger liquid [4] using the expression (8) for 
Xk- In that case, Eq. (20) reads 

rPu> D /2 J( x ) 

1 = VN(0) / dx nr l f \ T/ r , (21) 

V J Jo l-(3T{p)I{xY v ; 

where T(p) = [af(p) + cos ^)}(t') 2 /t, and 

tanh(xu) 



, « 1 f 1 t&niiixu) 
J l x ) = — / du 



In the limit f3u)D "C 1 the last integral is of elliptic form, which has a property 
fa dzj (y/z — a\Jb — z) = 7r; hence I(x) = 1. In effect, we obtain a remarkably simple 
expression for T c 



T c = WV(0)^ + i (t 



2 2 t 



T c » + T' c cos (-^j . (22) 



This result is of the same type as the formula derived for a = in the Fermi liquid case 
[1] (it can be extended to the case a ^ in straightforward manner). Hence, the transition 
temperature expression is universal and contains a combined effect of the intrinsic in-plane 
pairing and the two interplanar Josephson tunneling process.es 

To test the relative role of the two tunneling contributions we have fitted the tabulated 
[8] T c values to the expression (22). Table I displays the detailed comparison, together with 
the T c ° and T' c values. The agreement between the theory and experiment is excellent for 
p > 1 systems. For p — 1 systems we have chosen T c ° as the middle point of the allowed 
interval; at least part of the T c ° value must be ascribed to the interplanar tunneling ~ 
In connection with this, one can say that to have a proper behavior in the limit of large 
p one must assume that T(k) ~ p^ 1 ; then af(p) ~ 1/p, and the contribution due to the 
intergroup tunneling is substantially smaller for larger p. This is the reason why T c ° value is 
essentially constant for p > 1. From this point of view, uncertainty in the T c values for the 
p — 1 systems can be understood by the planes intersection due to the sample imperfection 
thus hampering the tunneling. Also, because of T c ° independence on the systems with larger 
p, the coherence length along the axis must exceed the smallest distance between the two 
neighboring groups. Finally, the maximal available T c for each family is for p — > oo, and is 
also listed in Table I. 

From our results, one can draw a very important prediction: the T c should increase 
substantially by applying the pressure along the c axis, as the pressure will increase it 
towards the value T c (oo). 

In summary, we have derived analytically the explicit forms of the space profile of the 
superconducting gap, and the transition temperature; these results were obtained without 
the need of specifying the pairing potential in the plane, and independently of the nature of 
the normal state. 
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TABLES 



TABLE I. Critial temperatures for various single and multilayer superconductors 



Material 


P 


rpexp 


7^0 

c 


T' c 


Theory: T c 


Bi 2 Sr 2 Cu0 6 


1 


0-20 


10 


148 


10 


Bi 2 Sr 2 CaCu0 8 


2 


85 


10 


148 


84 


Bi 2 Sr 2 Ca 2 Cu 3 Oio 


3 


110 


10 


148 


112 


— 


oo 


- 


10 


148 


158 


Tl 2 Ba 2 Cu0 6 


1 


0-80 


41 


126 


41 


Tl 2 Ba 2 CaCu0 8 


2 


108 


41 


126 


103 


Tl 2 Ba 2 Ca 2 Cu 3 Oi 


3 


125 


41 


126 


128 


— 


oo 


- 


41 


126 


167 


Tl 2 Ba 2 Cu0 5 


1 


0-50 


24 


121 


24 


Tl 2 Ba 2 CaCu 2 7 


2 


80 


24 


121 


84 


Tl 2 Ba 2 Ca 2 Cu 3 9 


3 


110 


24 


121 


107 


Tl 2 Ba 2 Ca 3 Cu 4 0ii 


4 


122 


24 


121 


122 




oo 




24 


121 


145 


HgBa 2 Cu0 4 


1 


94 


94 


59 


94 


HgBa 2 Ca 2 cu 3 8 


3 


135 


94 


59 


135 




oo 




94 


59 


153 
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FIGURES 

Fig. 1. Schematic plot of the gap magnitude (top, a) for a multilayer situation (bottom, 
b) involving groups of p > 1 tightly spaced planes. 
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